Resonant emergence of global and local spatiotemporal order in a nonlinear field 

model 
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We investigate the nonequilibrium evolution of a scalar field in (2+1) dimensions. The field is set 
in a double-well potential in contact (open) or not (closed) with a heat bath. For closed systems, 
we observe the synchronized emergence of coherent spatiotemporal configurations, identified with 
oscillons. This initial global ordering degenerates into localized order until all oscillons disappear. 
We show that the synchronization is driven by resonant parametric oscillations of the field's zero 
mode and that local ordering is only possible outside equipartition. None of these orderings occur 
for open systems. 
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The emergence of spatiotemporal ordered structures 
in nonlinear systems is an ideal laboratory for investigat- 
ing the trend toward complexification observed in nature 
at the physical, chemical, and biological level 0]. For 
these ordered structures to survive, they must interact 
with an external environment, which maintains the lo- 
cal nonequilibrium conditions. Examples can be found 
in hydrodynamics, in networks of chemical reactions 0, 
and in living organisms 3]. 

In field theory and cosmology, most of the interest in 
ordered configurations has focused on topological or non- 
topological static solutions of the equations of motion 4] . 
An exception to this trend are oscillons, long-lived time- 
dependent localized field configurations that have been 
found in field theory soft condensed-matter sys- 

tems j2j, and stellar interiors 0- Here, we show how 
oscillons spontaneously emerge as a nonlinear scalar field 
approaches thermal equilibrium. We observe not only the 
local emergence of spatiotemporal order (oscillons), but 
also that this emergence is initially synchronized (global) . 
Our results are applicable to any system modeled by a 
scalar order parameter cj) with amplitude-dependent non- 
linearities, if its potential V{(l)) satisfies d^V/d(p^ < for 
a range of (/). 

Consider a (2+1) dimensional real scalar field (j){x,t), 
which evolves under the equation of motion 



dw 



(1) 



The double-well potential Vdw(0) = -1(0^ — 1)^ has min- 
ima at = ±1. It has been shown that this system can 
generate oscillons, characterized by a persistent oscilla- 
tory behavior at their core To see this, prepare the 
field with a Gaussian profile, (j){r) = (pa exp(— r^/i?^) — 1, 
and let it evolve via Eq. Q). All that is needed is that 
V^^{(j)a) < and that the initial radius R be larger than 
a bifurcation value i?osc [Ml^- The properties of these so- 
called deterministic oscillons have been extensively stud- 
ied in two j^, [l3| and three dimensions 0, Q| . 

We investigate open (canonical) and closed (micro- 



canonical) systems. In both cases, the field is initially 
thermafized in a single-well potential Vsw{(l>) — {(j) + 1)^, 
symmetric about (/) — —1, chosen so that Vj,.j,,(— 1) = 
V^^{—1). The thermalization is achieved by coupling the 
field to an external heat bath via a Langevin equation 



d^ct) d(t) . , 



-Kw('^)+C, 



(2) 



where the viscosity coefficient 7 is related to the 
stochastic force of zero mean ^(x, t) by the fluctuation- 
dissipation relation (fc^ = 1 and T is the temperature of 
the heat bath), (^(x, t)^(x', <')) = 2-fTS^{x-x')S{t -t'). 
The numerical evolution is implemented on a square lat- 
tice with periodic boundary conditions, using a stag- 
gered leapfrog method. Integrations are performed with 
Sx — 0.1 and St = 0.01 with 1024 lattice sites per side. 
The coupling to the heat bath continues until equiparti- 
tion is satisfied at temperature T. Results are ensemble 
averages from 50 independent realizations, which, given 
the large number of degrees of freedom, have energy E 
within AE/E < 10"^. 

The potential is then switched from Vsw to Vdw, and 
the system is tossed again out of equilibrium. This switch 
can be interpreted as an instantaneous quench on the 
system, implemented by varying some control param- 
eter a, as is customarily done in many applications. 
To see this, shift (f> — > + 1 so that Vsw — 4''^ and 
Vdw = l/4[02(a0 - 2)^\a=i. Note that when a = 0, 
Vdw = Vsw The quench then occurs instantaneously by 
sending a = to a = 1 (or equivalently, Vsw Kiw)- 

For the closed system the coupling to the bath is re- 
moved at the time of the switch (by setting 7 = 0), while 
for the open system it is kept on throughout the simula- 
tion. After the switch occurs, energy exchange between 
the nonlinearly coupled modes will again drive the system 
to equipartition. (We note that the initial and final equi- 
librium temperatures never differ more than ~ 0.15%. 
Thus, each experiment is referred to by its initial ther- 
malization temperature.) For open systems, equilibrium 
is attained within a time scale of 0(7""'^), without any 
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FIG. 1: Bottom: lattice-averaged field for various tempera- 
tures. Top: fraction of the field above the inflection point. 

emergence of ordered configurations. This is in marked 
contrast with closed systems. 

Figure H (lower half) shows the area-averaged field 
0av (t) for closed systems at various initialization temper- 
atures. The time is set so that the switch to the double- 
well potential occurs at i = 0. Note that (j)av{t) displays 
damped oscillations, indicating that the zero mode of the 
field transfers its energy to higher modes until it reaches 
its temperature-dependent final equilibrium value {4>)t- 
At low temperatures, such as T = 0.1, the oscillations 
have small amplitude and remain nearly periodic. At 
higher temperatures, such as T = 0.2, 4>av{t) oscillates 
beyond the left inflection point and fluctuations in the 
field probe the unstable region of the double-well poten- 
tial. This can also be seen in the upper half of Fig. ^ 
where the time evolution of the fraction of the field above 
the inflection point fin{{t) is shown. For temperatures 
above T ~ 0.185, over half the field probes the unsta- 
ble region. For even larger temperatures, T > 0.25, the 
whole field goes above the inflection point, signaling the 
approach to criticality. Above this temperature the field 
separates into large, slowly evolving thin- walled domains. 
(The critical temperature of this system, where (/)av — > 0, 
is Tc ^ 0.270 ±0.005.) 

Figure |21 shows a sequence of snapshots of the field at 
T = 0.2, spanning in time about half an oscillation pe- 
riod of an oscillon. In order to generate this figure and 
relate the emergent configurations to oscillons, the field 
is smoothed with an optimal (Wiener) filter • This fil- 
tering technique is especially useful in our situation, since 
its implementation requires knowledge of the unwanted 
thermal noise in the system. With our choice of thermal 
initial conditions, Boltzmann statistics provides us with 



FIG. 2: Snapshots displaying synchronous emergence of os- 
cillons in the 2-dimensional field, starting at t = 15.5 and 
spanning half the period of an oscillon. Simulations can be 
viewed in Ref. [l^ . 

the power spectrum at the beginning of the experiment, 
(|(/)(k, t = 0)p) = r/(fc2 + 2), which is to be removed by 
the filter. Throughout the experiment, only modes with 
< |k| < 0.8 amplify above the noise, sometimes by as 
much as two orders of magnitude. Their spectrum has 
a shape and width that coincide with those of determin- 
istic oscillons. The transformation back into real space 
clearly reveals the emergence of localized field configura- 
tions within a smooth background. 

With the filtered field, we can catalog and track the 
location of all local extrema at each instant in time, the 
vast majority of which are seen to correspond to the cen- 
ters of localized configurations. We then measure the 
value of the field at each extrema and the correspond- 
ing configuration radius (at half maximum). With am- 
ple sorting, we compile a library containing each large- 
amplitude fluctuation during the entire evolution of the 
system {0 < t < 1500), from which we can obtain their 
nucleation times, sizes, periods of oscillation, and life- 
times. A deterministic oscillon is characterized both 
by its large-amplitude oscillations (the field at its core 
probes the positive half of the potential, (pa > 1) and its 
extreme longevity. We thus establish two criteria to se- 
lect the subset of all configurations which correspond to 
oscillons: they must have a maximum amplitude above 
the background satisfying (jja > ^ and they must survive 
for at least ten oscillations {t > 60). 

The inset in Fig. 13 shows the probability distribution 
functions for the radii R and periods of oscillation P for 
all the oscillons present throughout the simulations, with 
binwidth 5R — 6P = 0.1. The fitted curves are Gaussian 
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FIG. 3: The number of oscillons nucleated between t and t+5t 
at T — 0.2, with St = 1. The global emergence is evident early 
in the simulations. Inset: the probability distribution of radii 
and periods of oscillation. 



FIG. 4: Lines of constant amplification rate for small- 
amplitude modes at various temperatures, beginning with 
77^^ ~ 27 for the bottom-most contour and increasing in in- 
crements of Jt? = 0.05. 



functions with centers _Ro = 2.86 and Pq — 6.12 and 
widths CT_R = 0.33 and ap = 1.19, respectively. 

In Fig. O we also show the distribution of nucleation 
times for these oscillons at T = 0.2. This function gives 
the number of oscillons nucleated between t and t + St, 
with St = 1. The sharp peaks at early times, t < 60, cor- 
respond to the synchronous emergence of oscillons, while 
for t > 60 this global ordering gives way to only local 
ordering, in which oscillons emerge at arbitrary times 
with similar probability. Even this local emergence dis- 
appears after approximately t > 500 (discussed below). 
Notice the correlation between the nucleation activity of 
oscillons and the energy loss from the zero mode (Fig. 
^ at early times. This is observed for all temperatures 
within 0.185 < T < 0.25. 

To understand the origin of this synchronous emer- 
gence, we decompose the field as <^(x, i) — (/'av(0 + 
S(j>{x, t) and first investigate the behavior of 4>a.-v{t): prior 
to relinquishing its energy to higher modes, by approx- 
imating S(j){x., t) as remaining statistically equivalent to 
its initial symmetric state. Upon substituting this form 
of 0(x, t) into Eq. (Q) and performing a statistical av- 
erage over the fluctuations, we arrive at the mean-field 
equation of motion for (/>av (i) : 0av = [1 — 3(^(/)^)]0av — 0av' 
where (50^) is the two-point correlation function at t = 
and depends linearly on the temperature T. Thus, <piiv{t) 
starts at = — 1 and oscillates anharmonically in the 
left-hand well of an effective double-well potential, with 
minima at (j)± = ±-\/l — 3{S(f)'^). This mean-field ap- 
proximation works very well at describing the evolution 
of 4'avit) at temperatures sufficiently far from Tc (c.f. 



T = 0.1 in Fig. It also gives (0)t « to within 4% 
even at the highest temperatures we considered. 

We now examine the behavior of small fiuctuations 
about (/)av (t) • Linearizing Eq. l(T]l with respect to S6(x. t) 
and taking the Fourier transform we obtain (for fc > 0) 

5>(A:, t) + + Kiw (0av(t))] S(t>{k, t) = . (3) 

Equations of this type, generalized Mathieu equations, 
are known to exhibit parametric resonance, which can 
lead to exponential amplification expryt) in the os- 
cillations of S(j){k,t) at certain wavelengths, in response 
to the time-dependent harmonic term. These equations 
have been of great interest in reheating studies of infia- 
tionary cosmologies . To verify that this is the mech- 
anism behind the synchronous amplification of oscillon 
modes early in the simulations, we use in Eq. (j^l the an- 
harmonic solution for 0av(O obtained in the mean-field 
approximation above. 

Figure ^ shows lines of constant amplification rate rj 
of the fluctuations S<t){k, t) for various k and T. At 
low temperatures, T < 0.165, no modes are ever am- 
plified. As the temperature is increased, so is the am- 
plitude of oscillation in 0av(i), eventually causing the 
band 0.6 ^ k < 0.9 to resonate. We note that the 
characteristic wavenumber of oscillon configurations is 
A:osc = 2/i?o Using the result in Fig. Hfor the aver- 
age configuration radius Rq — 2.86, we obtain /cqsc — 0.7. 
From Fig. 0] modes with kosc — 0.7 are excited for 
T > 0.18, the temperatures above which we see the 
synchronized emergence of oscillons. For T > 0.22, a 
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FIG. 5: The change of Tl{t) from the initial state for closed 
systems at T — 0.2. The exponential approach to equilibrium 
is clear at late times. The inset iUustrates the role of oscillons 
as a bottleneck to equipartition. 



n(t) (dotted line) that arise due to the synchronous os- 
cillations in the kinetic energy of these configurations. 
Also shown (solid line) is the average between successive 
peaks of n(t), with a plateau at 20 < t < 70 that coin- 
cides with the maximum oscillon presence in the system. 
Thus, oscillon configurations serve as early bottlenecks 
to equipartition, temporarily suppressing the diffusion of 
energy from low (0 < |k| < 0.8) to higher modes. 

We have investigated the nonequilibrium evolution of 
a scalar field with a double- well potential. For a range of 
temperatures, the approach to equilibrium is character- 
ized by three stages: first, the synchronized emergence of 
oscillons; second, the loss of the initial synchronicity, but 
the persistence of oscillons; and third, their disappear- 
ance as the system approaches equipartition. It would 
be interesting to investigate the possibility of controlling 
the duration of the synchronization stage and search for 
this emergent behavior in laboratory systems, ranging 
from vibrating grains to ferromagnetic materials. 

MG was supported in part by NSF grants PHY- 
0070554 and PHY-0099543. 



second band of longer wavelength modes becomes ex- 
cited as well, signaling the onset of criticality. These 
general results corroborate our earlier findings: for tem- 
peratures 0.185 ^ T < 0.25, oscillations of the zero 
mode drive, via parametric resonance, amplification of 
the modes comprising oscillon configurations, ultimately 
leading to the usual breakdown of the linear approxima- 
tion. The growth of these modes corresponds to large- 
amplitude fluctuations that probe the unstable regions of 
the potential {d'^V/d4>^ < 0), which coordinate to form 
the observed coherent structures. 

Finally, we introduce a measure of the partitioning of 
the kinetic energy n(f), which we use to describe the 
nonequilibrium evolution of the system: 



n(i) ^ ~ J (fk p(k, t) lnp(k, t), 



(4) 



where p(k,t) = K{k,t)/ J (fkK(k,t) and K{k,t) is the 
kinetic energy of the fcth mode. Hit) attains its maxi- 
mum [Ilinax ~ In(iV) on a lattice with TV degrees of free- 
dom] when equipartition is satisfied. This occurs both 
at the initial thermalization (t = 0) and final equilib- 
rium states, since in this case all modes carry the same 
fractional kinetic energy. In Fig. Elwe show the change 
of n(t) from the initial state, n(i = 0) — n(t), for the 
closed system at T = 0.2. At late times {t > 150), we 
have found that the system equilibrates exponentially in 
a time scale r ~ 500. At early times, the localization 
of energy at lower k modes, corresponding to the global 
emergence of oscillons, prolongs this approach to equipar- 
tition. The inset of Fig. [S] shows the large variations in 
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